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Nurse  reus  articles  have  been  written  concerning  the  rorticity  equation  and 
the  computation  of  constant  absolute  vorticity  trajectories.  The  purpose  of 
unis  investigation  is  to  determine  statistically  the  degree  of  accuracy  of 
these  computed  trajectories  for  twenty-four  and  forty-eight  hour  periods,  and 
if  possible,  to  find  some  systematic  departure  of  the  trajectories  so  that  a 
correction  factor  might  be  applied  to  computations  made  in  the  future 


Methods 

The  Bellamy  trajectory  slide  rule  was  not  used  for  making  the  c imputations 
tested  because  it  is  based  on  the  assumption  of  the  earth  as  a plane.  Instead, 
tables  were  used  which  were  prepared  by  interpolation  of  data  taken  from  curves 
made  by  the  differential  analyzer  at  the  WBAH  Analysis  Center.  These  assume  a 
spherical  earth.-1 

Several  points  for  computation  were  selected  each  day  from  the  500-«b 
chart,  with  the  criteria  for  their  selection  being  the  same  *a  established  in 
the  theory,  i.e.,  at  inflection  points  or  in  straight-line  flow,  in  art>as  of 
no  shear  in  the  horizontal  and  in  areas  of  no  acceleration  along  the  stream- 
lines. 

The  wind  direction  and  speed  were  taken  from  actual  data  whenever  possible 
or  were  computed  using  the  geostrophic  assumption.  Since  non-divergence  is 
another  assumption  of  the  conservation  of  vorticity  theorem,  tho  vied  was 
reduced  to  the  level  generally  considered  to  be  the  level  of  con-divergence, 
i.e.  600  mb.  The  method  of  reduction  was  to  subtract  two-thiids  of  the  vertical 
wind  shear  between  700  mb  and  !>00  r«b  from  the  speed  obtained  at  £00  mb.  The 
wind  direction  was  not  changed.  The  tables  *'sre  then  entered  with  these  values 
for  the  wind  and  the  four  parameters  descrioing  the  trajectory  vere  obtained. 


The  position,  wind  direction  and  velocity,  and  the  height  of  the  500~nb 
surface  at  the  point  selected  were  recorded,  as  was  the  prognostic  and  veri- 
fication values  of  wind  direction,  and  the  height  of  the  500-nb  surface.  The 
verification  wind  direction  was  taken  from  the  500-rab  contours  at  the  point  in 
question.  Separate  records  were  kept  of  twenty-four  and  forty-eight  hour 
verifications.  Wind  speed  was  not  verified.  In  all,  seventy-five  trajectories 
were  computed  during  the  period  mid-October  to  mid-December  19>1. 

1.  The  personnel  involved  were  changed  during  the  test  and  were  of 
widely  varying  experience  in  the  making  of  these  computations. 

2.  A standard  reduction  of  wind  speed  from  SCO  mb  to  600  mb  roust,  at 
times,  give  erroneous  values  of  wind  speed..  . 


The  tables  and  instructions  are  appended. 
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3.  The  points  were  selected  from  a hemispheric  chart  with  a polar  pro- 
jection. This  projection  makes  the  selection  of  the  inflection  point  extremely 
difficult. 

U.  Since  the  trajectories  were  primarily  for  use  in  a laboratoiy  course, 
the  point  of  computation  was  selected  su*5*  that  the  verifying  wind  would  be 
in  or  near  the  forecast  area  in  forty-e4  hours.  This  necessitated  the 
initial  point  of  many  of  the  trajectories  t.  be  taken  in  the  East  Pacific  area 
where  analysis  is  at  best  not  precise. 

Examination  of  results 


Table  1 gives  the  departure,  in  10 ’s  of  degrees,  of  the  observed  wind 
direction  from  the  forecast  direction^  with  positive  values  indicating  an 
observed  wind  to  the  right  (clockwise)  from  the  forecast  wind.  The  values  for 
the  2U-hour  period  show  a significantly  good  verification  with  a reasonably 
normal  distribution  and  about  half  of  the  cases  verifying  within  the  accuracy 
of  the  wind  measurement  (t  10°  from  aero).  The  values  for  1*8  hours,  on  the 
other  hand,  do  not  show  a normal  distribution  and  do  not  show  a pronounced 
maximum  near  the  zero  departure.  Titus  the  general  conclusion  that  such  tra- 
jectories have  considerable  value  for  a 2 It  -hour  period  but  little  value  for  a 
i»8~heur  period. 

Table  2 is  the  result  of  the  attempt  to  find  some  systematic  departure 
of  the  trajectories.  Only  ? labour  verifications  were  used  due  to  the  poor 
distribution  of  the  U8-hour  verifications.  The  various  groupings  were  selected 
before  examination  of  the  data. 

The  breakdown  of  initial  direction  was  roughly  winds  in  the  southwest 
quadrant  versus  those  in  the  northwest  quadrant,  or,  one  might  say,  winds 
approaching  a ridge  versus  winds  approaching  a trough.  These  distributions 
show  no  appreciable  difference  as  both  approximate  a normal  distribution. 
Another  choice  of  longitude  intervals  would  not  change  the  result  significantly 
as  the  points  were  quite  scattered. 

The  two  wind  speed  groupings  are  quite  similar  and,  therefore,  they  too 
cannot  be  used  to  obtain  a systematic  departure. 

The  groupings  for  the  initial  longitude,  or  the  other  hand,  do  show  some 
difference.  They  were  chosen  so  as  to  approximately  separate  continental 
from  oceanic  areas.  This  breakdown  can  be  classified  generally  as  the  data 
versus  no  data  excepting  the  region  near  the  stationary  ship  in  the  Gulf  of 
Alaska.  As  this  exception  would  affect  the  result  only  by  improving  the  "no 
data"  distribution,  it  can  readily  be  said  that  verification  will  be  better 
when  computations  are  made  either  using  actual  data  or  having  data  nearby. 

Not  a startling  fact,  but  one  worth  considering. 

The  verifying  longitude  also  shows  an  interesting  difference  in  the  three 
groups  selected.  Assuming  the  potential  vorticity  to  be  conserved,  winds 
should  deviate  to  the  right  of  the  CAVT  when  the  vertical  column  is  becoming 
less  and  vice  versa.  The  grouping  of  from  12$°W  to  110°W  was  selected  as  that 
area  where  the  air  column  would  have  shrunk  due  to  passage  over  the  mountains. 
The  105  W°to  100°W  group  is  Just  east  of  the  mountains  where  a minimum  devia- 
tion would  be  expected  and  the  95°W  to  80°W  group  is  where  the  column  should 
have  stretched  enough  to  show  deviation  to  the  left.  The  figures  hear  out 


-3- 


this  selection  to  sane  degree  although  the  flO°  deviatioi*  still  predominates. 

The  few  computations  at  higher  latitudes  were  removed  from  this  sample. 

The  final  classification  in  Table  2 was  selected  a3  it  was  felt  tnat  the 
change  in  the  height  of  the  particle  during  the  forecast  period  might  indicate 
energy  changes  of  tne  particle  and  thereby  give  a systematic  deviation.  This 
hypothesis  again  is  borne  out  to  some  extent  by  the  verification  as  the 
interval  1200  ft  gives  a much  more  normal  distribution  than  the  other  grouping.. 


A final  attempt  to  obtain  a prognostieally  valuable  systematic  deviation, 
or  raiher  non-deviation,  was  to  try  and  find  a group  of  successive  days  in 
which  the  verification  was  fairly  good.  The  data  cards  were  checked  for  a 
string  of  -20°  verifications  and  the  dates  recorded.  Three  groupings  of 
greater  than  two  days  shewed  up.  Ihey  were  as  follows? 


1.  22  October  through  1 November-17  trajectories  3 exceptions 

2.  6 November  through  12  November,  10  trajectories  1 exception 

3.  22  November  through  2 December,  12  trajectories  1 exception 

Tne  exceptions  are  verifications  slightly  larger  than  *20°.  These  dates  were 
bracketed  on  a time  section  of  the  aonal  profile,  as  confuted  at  Andrews  Field 
and  sent  out  on  the  circuit  n0"  teletype.  Three  day  overlapping  totals  were 
used  to  smooth  out  irregularities  in  the  daily  profile.  It  was  then  discovered 
that  each  of  the  brackets  was  centered  on  an  Nil  stage  of  the  index  cycle  (1). 
This  supports  the  notion  that  vorticity  is  conserved  to  a great  extend  during 
high  index  for  this  index  always  occurs  diring  an  Nil  stage. 


Summary 


In  sp ite  of  the  factors  noted  which  might  detract  from  the  results,  the 
following  was  shown  more  or  less  significantly  concerning  constant  vorticity 
trajectories  computed  under  the  assumption  of  a spherical  earth* 

1.  They  do  not  verify  well  after  a U8-hour  period. 

2.  They  will  verify  well  after  a 2h-hour  period. 

a.  They  worked  as  well  an  yv  fie  re  within  the  ranee  of  speeds 
and  direction  tested  (30  to  80  knots  and  170°’  to  360°). 

b.  The  verification  was  somewhat  betler  for  starting  points 
over  the  North  American  continent  than  for  points  over  the 
ocean  (probably  due  to  better  data  and  malysis). 

c.  The  systematic  deviation  tends  to  be  to  the  right  over 
maintain  areas  and  io  the  left  soma  distance  to  the  lee. 

d.  Successive  days  with  good  verification  tend  to  occur 
during  the  Nl'I  stage  of  the  index  cycle. 


(1)  Riehl  et  al,  l?f>2.  Forecasting  in  Middle  Latitudes,  A.M.S.  Monograph 
No.  5* 


CAVT  Tables 


The  data  contained  in  these  tables  were  extracted  graphically  from  CAVT 
trajectories  made  by  the  differential  analyser  at  the  HBAH  Analysis  Center* 

The  trajectories  used  were  for  latitudes  40°,  50°,  60°  and  70°N,  wind  direotions 
of  250°,  210°  tnJ  170°,  and  wind  speeds  of  30,  50  and  70  knots , making  36 
trajectories  in  all*  To  obtain  the  tables , a linear  interpolation  was  assumed 
to  be  valid*  For  this  reason,  those  tables  do  not  pretend  to  be  striotly 
oorrect*  They  are,  however,  a refinement  of  the  Bellamy  Slide  Rule  trajeotories , 
simc  the  earth’s  curvature  is  considered  in  this  method* 

The  tables  require  the  same  parameters  as  the  Bellamy  Slide  Rule)  that 
is,  latitude,  wind  direction,  and  velocity*  The  values  then  obtained  art  as 
follows  f 

Dy  - Distance  in  degrees  loigitude  between  the  inflection  points  of 
the  upper  loop* 

Yy  ■ Amplitude  in  degrees  latitude  of  the  upper  loop* 

m Distance  in  degrees  longitude  between  the  inflection  points  of 
the  lower  loop* 

* Amplitude  in  degrees  latitude  of  the  lower  loop* 

To  obtain  the  position  of  the  particle  at  some  future  time,  it  is  neoes- 
s ary  to  multiply  the  wind  speed  by  the  time  interval  desired  and  lay  this 
dlstanoe  out  along  the  trajectory* 

Notes  For  wind  direotions  250°  through  170°  oounterolookwise,  draw  arrow 
through  origin  and  measure  angle  counterclockwise  from  90°,  e*g*  for  wind 
direction  230°  enter  table  using  OC-  40°,  For  wind  direotions  290°  through 
010°  olookwise,  draw  arrow  through  origin  and  measure  angle  clockwise  from 
90°*  Then  plot  lower  loop  of  trajectory  first,  e,g*  for  wind  direction  310° 
enter  the  table  using^  - 40°, 
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Table  1— Departure  of  Observed  Wind  Direction  fro* 

Forecast  Wind  Direction  (as  per  cent  of  total) 


Table  2 — 2u  hour  Departure  of  Observed  from  Forecast  Wind 
Direction  vs.  Various  Paraaeters  (as  per  cent  of 
totJ.) 


